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Abstract 

The production rate for top quarks at the Fermilab Tevatron is presented 
using the exact order corrected cross section and the resummation of the 
leading soft gluon corrections in all orders of perturbation theory. 
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1 Introduction 



At the Tevatron, the top quark will be mainly produced through ti pair 
creation. Both top quarks will then decay to (W, b) pairs, and then each W 
can decay either hadronically or leptonically. It is in the channel where both 
W's decay leptonically, one to a (e, v e ) pair, the other to a (//, z/ M ) pair, that 
a large part of the current search effort to find the top quark is concentrated. 
This is because the background in this channel from W + W~ plus jets is fairly 
small. Using the known branching fractions of the above decays, and taking 
acceptances into account, an experimental top production cross section is 
determined, which is then compared with a curve of the theoretical cross 
section as a function of the top quark mass m. At present the top quark 
has not yet been discovered at the Tevatron, and thus an estimate of the 
theoretical cross section is needed to either determine its mass or establish 
a lower limit on it. We therefore discuss the estimates of this cross section 
here. 

There have been previous predictions of the top quark cross section based 
on the results of the fixed order a 2 s plus otl contributions in perturbative QCD 
(pQCD) 0,@, @,g. As with all fixed order pQCD calculations, these 
contain a scale (factorization scale = renormalization scale) which reflects 
the size of the uncalculated 0(af) and higher order terms. Although the 
dependence on this scale is relatively flat, indicating that the present result 
is stable under scale changes, the size of the O(af) term is disturbing. In 
fig.l we show the ratio R = (cx^Lt + CT exLt)/ CT cxact f° r top quark masses 
in the range relevant to the Fermilab Tevatron. With a^ ct we denote the 
exact order a^ n+2 ^ contribution to the cross section, which implies that a^ct 
stands for the Born contribution. Here, and throughout this paper, we have 
used recent parton distribution functions (MRSDJ) | ||, choosing the DIS 
factorization scheme, and the two-loop running coupling constant (in the 
MS scheme) with five active flavours and A = 0.152 GeV. The ratio R of 
the next-to-leading order (NLO) to the leading-order (LO) term is usually 
referred to as the K-factor. 

One notes from fig.l that the higher order corrections in the qq channel 
are small, whereas those in the gg channel are 70% or larger. Throughout the 

2 We thank W.J. Stirling for sending us the DIS scheme MRS parton distribution 
fruitions. 
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mass range we consider the qg and qg channel give negligible contributions so 
we do not consider them. To show which channel has the largest cross section 
in this mass range we plot the ratios of the NLO gg and qq contributions to 
the total result in fig. 2. Here we see that as the top quark mass increases 
the qq channel contribution is larger than the gg channel one, so the effect of 
the large K-factor in the latter channel, as seen in fig. 1, decreases. However, 
at masses around 150 GeV/c 2 the gg channel still contributes around 20%. 
This is enough to worry about even higher order corrections. 

These large corrections are predominantly from the threshold region for 
heavy quark production, where we have shown previously that initial state 
gluon bremsstrahlung (ISGB) is responsible for the large corrections at NLO 
JJ. To demonstrate this we have run our cross sections programs with a cut 
on the variable rj = (s— 4m 2 ) /4m 2 , where s is the square of the parton-parton 
cms energy so that: 

fVcut I 

v(Vcut) = H / dr lT— $ y(^>M 2 ) Vij(v,™ 2 ^ 2 )> (1) 

where $y denotes the parton flux and &ij the partonic cross section for the 
incoming partons % and j. The precise definitions of these functions are 
given in Ref. ||. The factorization scale is fi, which we have set equal to the 
renormalization scale. Notice that the maximum value of r] cut is given by 
Vmax = (S — 4m 2 ) /4m 2 , where S denotes the square of the total hadronic 
energy in the cms system. As we increase r) cut from small values (~ 10 _1 ) to 
larger values, where the actual value of the cross section is approached, there 
is a rapid rise in a(f] cut ). Figures 3 and 4 show cr(r] cut ) for the gg and the 
qq contributions to the cross section. The fact that both cross sections rise 
sharply around 7] cut = 1 (where s ~ 8m 2 ) indicates that the threshold region 
is very important. This is especially true in the gg channel, which dominates 
the total cross section at smaller top quark masses. Both figures also show 
that the cross section flattens out if r] cut is increased further, indicating that 
partonic processes with s >> 4m 2 contribute very little to the cross section. 

In a previous paper j?| we carefully examined the dominant logarithms 
from ISGB which are the cause of the large corrections near threshold. Such 
logarithms have been studied previously in Drell-Yan (DY) |§ production at 
fixed target energies (again near threshold) where they are responsible for 
correspondingly large corrections. In [|7| we exploited the analogy between 
DY and heavy quark production cross sections and proposed a formula to 
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resum the leading and next-to-leading logarithms in pQCD to all orders. 
Since the contributions due to these logarithms are positive (when \x = m), 
the effect of summing the higher order corrections increases the top quark 
production cross section over that predicted in O(a^). This sum, which will 
be indentified with <7 res , depends on a nonperturbative parameter /x . The 
reason that a new parameter has to be introduced is due to the fact that the 
resummation is sensitive to the scale at which pQCD breaks down. As we 
approach the threshold region other, nonperturbative, physics plays a role 
(higher twist, bound states, etc) indicated by a dramatic increase in a s . We 
chose to simply stop the resummation at a specific scale /i where A << 
fiQ << m since it is not obvious how to incorporate the nonperturbative 
physics. Note that our resummed corrections diverge for small /xo but this is 
not physical since they should be joined smoothly onto some nonperturbative 
prescription and the total cross section will be finite. However, at the moment 
our total resummed corrections depend on the parameter fi Q for which we can 
only make a rough estimate. See for more details. 

Let us begin by showing the effects of including only the leading soft gluon 
contribution at 0(af), which we call a® . An explicit expression for this 

contribution is given in Here a^ p stands for the approximation to c^Lt 
where only the leading soft gluon corrections are taken into account. Figure 
5 shows three curves for the exact cross section calculated through 0(a%) at 
the scales \i = 2m, m and m/2. This is the traditional method of estimating 
the size of uncalculated higher order contributions. For comparison we add 
to the ii — m case the approximate 0(atf) contribution, yielding a total cross 
section which is not in the range spanned by the previous curves, but slightly 
above. Therefore, the traditional method does not work very well in this case, 
due to the size of the corrections, as was already pointed out in || . 

Now we study the effect of the resummation, which depends on /i , by 
calculating <r res . However, because we know the exact 0(a^) result, we can 
make an even better estimate of the cross section by calculating the quantity 

<^imp <^res ^app ^cxact; (2) 

which we call the improved cross section. We remind the reader that cr*™-* 
denotes the 0(a^ n+2 ^) contribution to the cross section. Further a^ct denotes 
the exact calculated cross section and a™ the approximated one where only 
the leading soft gluon corrections are taken into account. We compare the 
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results for cr imp versus the fixed order result OexLt + ^exLt + CT ipp ^ n for 
various (reasonable) values of /io- Note that fiQ need not be the same in the 
qq and gg channels because the convergence properties of perturbation series 
could be different in these channels and depend on the factorization scheme. 
For example the cross section due to the qq process seems to converge faster. 

As we decrease /jlq the cross sections increase. In figure. 6 we show three 
curves for various choices of /jlq. The requirement that A << /x << m t is 
satisfied by all three choices. This is not a very restrictive requirement in the 
sense that it still leaves a large range of values possible, thus rendering the 
cross section from pQCD more uncertain than thought previously. 

In fig. 7 we show the contribution from the qq channel to fig. 6, while the 
gg channel is given in fig. 8. The plots show that the range of possible cross 
sections is quite narrow for the qq channel (due to the smaller relative correc- 
tion) while it is relatively large in the gg channel. The curves are calculated 
in the DIS factorization scheme, where the corrections are smaller because 
most of them have been absorbed in the parton distribution functions. Thus 
the resummation is successful for the qq channel, in the sense that o"i mp dif- 
fers very little from cx^Lt + CexLt- This is unfortunately not the case for gg. 
We have also checked that these results change very little when using CTEQ 
rpTOf] parton distribution functions. This is to be expected because, due to 
the fact that top is so heavy, the cross section is mainly sensitive to parton 
distribution functions at large x where they have been well measured. 

Finally in Table 1 we present a lower limit estimate (from |[L1||), and a 
central value and upper limit estimate. The latter two are also shown in fig. 6 
as the central and upper solid line respectively, and are obtained as follows. 
For both estimates we used the improved cross section (|2]). Each of the three 
terms in @ was calculated according to ([[]), with r\ cut = r\ max . Furthermore, 
as stated earlier, we used the MRSDJ parton distribution functions, and 
chose the DIS factorization scheme and the two-loop running coupling con- 
stant (in the MS scheme) with five active flavours and A = 0.152 GeV. The 
exact partonic cross sections used to calculate cx^Lt + 0"exLt were obtained 
from [p], ^ H|]. The approximate DIS scheme partonic cross section used for 
determining is given explicitly in eqn. (2.10) in |7]]. For the partonic re- 
summed cross section we used eqn. (3.24) in the same reference. For all three 
cross section contributions we chose // = m. However, the central value and 
upper limit estimates differ by the value chosen for the nonperturbative pa- 
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rameter /xq. Thus, for the central value we chose /xo = 0.1m and /xo = 0.25m 
for the qq and gg channels channels respectively, whereas for the upper limit 
we chose /x = 0.05m and /x = 0.2m, respectively. It should be noted that 
these estimates merely represent what we think are reasonable choices for 
the parameter /x , and are thus still not completely rigorous. As far as the 
lower limit estimate is concerned, note that we can write 
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"imp <-> exact ^ "exact ^ / , "app' \°J 

i=2 

Since crW p > for all i at /x = m |7|], the true total cross section is likely 
larger than 

Slower = C"exact ^cxact °~app> (^0 

so we are justified in using this value as a lower limit. The calculation of the 
three terms in (f|) goes analogously as described above for the central value 
and upper limit cases. It differs only in the fact that we do not have the 
nonperturbative parameter /xq here and that we used the conservative value 
A = 0.105 GeV in the expression for a s when we make this estimate for the 
lower limit. The explicit expression for cr^ p in the DIS scheme is given in 
eqn. (2.14) in 0. 

In conclusion we have demonstrated that the leading and next-to-leading 
logarithmic corrections to the top quark production cross section are large 
near threshold and have to be resummed to give a more precise estimate for 
the cross section. When the ISGB contributions are resummed to all orders 
in pQCD a new scale /x has to be introduced, which measures the sensitivity 
of the cross section to nonperturbative physics. The top quark production 
cross section at the Fermilab Tevatron is sensitive to this new scale, mainly 
via the contribution from the gluon-gluon fusion channel. 
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Table 1. The first and fifth column contain the top quark mass in GeV/c 2 . 
The columns denoted by 'Lower' show our lower limit estimate of the top 
quark cross section in picobarns, the columns denoted by 'Central' show our 
central value estimate, and the columns denoted by 'Upper' show our upper 
limit estimate. 



7 



References 

[1] P. Nason, S. Dawson and R.K. Ellis, Nucl. Phys. B303 (1988) 607; B327 
(1989) 49; E. B335 (1990) 260. 

[2] G. Altarelli, M. Diemoz, G. Martinelli, and P. Nason, Nucl. Phys. B308 
(1988) 724. 

[3] W. Beenakker, H. Kuijf, W.L. van Neerven and J. Smith, Phys. Rev. 
D40 (1989) 54; W. Beenakker, W.L. van Neerven, R. Meng, G.A. Schuler 
and J. Smith, Nucl.Phys. B351 (1991) 507. 

[4] R.K. Ellis, Phys. Lett. B259 (1991) 492. 

[5] A.D. Martin, R.G. Roberts, W.J. Stirling, Phys. Lett. B306 (1993) 145. 

[6] R. Meng, G.A. Schuler, J. Smith and W.L. van Neerven, Nucl. Phys. 
B339 (1990) 325. 

[7] E. Laenen, J. Smith and W.L. van Neerven, Nucl. Phys. B369 (1992) 
543. 

[8] G. Sterman, Nucl. Phys. B281 (1987) 310; D. Appel, P. Mackenzie and 
G. Sterman, Nucl. Phys. B309 (1988) 259; S. Catani and L. Trentadue, 
Nucl. Phys. B327 (1989) 323, B353 (1991) 183; H. Contopanagos and 
G. Sterman, Nucl. Phys. B400 (1993) 211. 

[9] F.A. Berends, J.B. Tausk and W.T. Giele, Phys. Rev. D47 (1993) 2746. 

[10] J. Botts, J.G. Morfm, J.F. Owens, J. Qiu, W.-K. Tung, H. Weerts, Phys. 
Lett. B304 (1993) 159. 

[11] E. Laenen, preprint FERMILAB-Pub-93/155-T, June 1993. 



8 



Figure Captions 



Fig. 1. The ratio R ('K-factor') for the NLO exact top quark cross section 
as a function of the top quark mass. Plotted are the K-factors for 
the total cross section (solid line), and individually for the qq channel 
(long-dashed line) and gg channel (short-dashed line). 

Fig. 2. Fraction of qq channel (long-dashed line) and gg channel (short- 
dashed line) contribution to total NLO cross section as function of the 
top quark mass. 

Fig. 3 Cross section as function of i] cut (see eq.(l)) for qq channel. We used 
the DIS scheme MRSD_' parton distribution functions, and m = 100 
GeV/c 2 . Plotted are cr(al) (solid line) and cr(a^) (dashed line). 

Fig. 4 Same as fig. 3 but now for the gg channel. 

Fig. 5 The NLO exact cross section as a function of the top quark mass for 
three choices of scale: /i = m/2 (upper solid line), fi = m (central solid 
line) and fi = 2m (lower solid line), and the NLO exact cross section 
plus the 0(a A s ) contribution at /i = m, (dashed line). 

Fig. 6 The O(a^) cross section at fi — m (dashed line) and cr imp (eq. (2)) 
for three choices of scale /zq, the two numbers per line corresponding to 
the qq and gg channels respectively: 0.05 m/0.2 m (upper solid line), 
0.1 m/0.25 m (central solid line), 0.2 m/0.3 m (lower solid line). 

Fig. 7 Range of cross sections for qq channel only. The two solid lines span 
the exact NLO cross section for the range m/2 < fi < 2m, and the two 
dashed lines a uap for the range 0.05m < /i < 0.2m. 

Fig. 8 Range of cross sections for gg channel only. The two solid lines span 
the exact NLO cross section for the range m/2 < fi < 2m, and the two 
dashed lines a inip for the range 0.2m < /x < 0.3m. 
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